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Abstract 

The resistivity of a one-dimensional Josephson junction array at zero temper- 
ature is calculated by estimating the nucleation rate of quantum phase slips. 
We choose a certain collective coordinate which describes the nucleation pro- 
cess and estimate the corresponding effective mass. In the strong coupling 
regime, where Josephson coupling energy exceeds the charging energy, we 
calculate the nucleation rate by means of WKB method. Our estimation is in 
good agreement with recent experimental data. The superconductor-insulator 

transition point is also discussed. 
74.50.+r, 73.40.Gk, 64.60.Qb, 64.60.My 



1 



Typeset using REVTfjX 



Quantum fluctuations have much influence on the properties of one-dimensional systems 
at low temperatures. For example, in extremely thin superconducting wires which can be 
regarded as nearly one-dimensional, the quantum fluctuations of the phase of the order 
parameter nucleate phase slips, which result in finite resistance at zero temperature [|TJ. 
Similar effects are expected for a one- dimensional Josephson junction (ID JJ) array. In a 
ID JJ array, the phase difference and charge difference across a junction are canonically 
conjugate variables, so the charging effect is the origin of the quantum fluctuation of phase. 
When the system size becomes small and the charging energy large enough, nucleation occur 
frequently and at some critical point, the quantum phase transition from a phase-coherent 
(superconducting) state to a charge-ordered (insulating) state occurs. Theoretically, this 
Superconductor-Insulator (S-I) transition in a ID JJ array can be understood as some ver- 
sion of Kosterlitz-Thouless transition in a (1+1) dimensional classical spin model, extra 
dimension being imaginary time 

Recently, such S-I transition was experimentally investigated in detail with ID JJ arrays 
having different numbers of superconducting grains [Q. In that work, the S-I transition was 
actually observed for each N near the point predicted by Kosterlitz-Thouless theory, and 
what is also remarkable, it was found that arrays have finite resistance at zero temperature 
even in the superconducting side. It is probable that these intrinsic resistance in the su- 
perconducting side are caused by the quantum nucleation of phase slips, as in the case of 
thin superconducting wires. In order to confirm this interpretation, we study the quantum 
nucleation process of a phase slip and calculate the nucleation rate in this Letter. The S-I 
transition point is also discussed. 

To begin with, we consider a linear array of iV(3> 1) superconducting grains embedded 
in an insulator at zero temperature. The grains are assumed to be small compared with 
bulk coherence length, so that the state of ith grain can be described by a single phase 
and number n« of Cooper pairs on the ith grain. The Euclidean action and Hamiltonian are 

SE[{niA}] = Jdr ^ihnJ^ + HiiniA})^, (1) 
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In Eq. (Q), we take into account two kinds of charging energies, Eq and Eq, the former 
is between a grain and the ground and the latter between nearest neighbor grains. Ej is 
Josephson coupling energy and n is the average number of Cooper pairs on a grain. Although 
the normal component is not explicitly mentioned here, the property of the system largely 
depends on the value of normal tunnel resistivity 0-I7J. In this Letter, we restrict our 
attention to the superconducting regime, so Rt is assumed to be smaller than the quantum 
resistance Rq> = ~ 6.45M7. 

A supercurrent state of the above Hamiltonian is metastable, even at zero temperature. 
Such a state corresponds to a local minimum of the energy functional of rii and 9{, and 
definitely decay to lower local minima, due to quantum tunneling. When the effective 
inertial mass of the phase slip is not greatly dependent on the tunneling path, the most 
probable tunneling path will be through the valley, passing through two local minima and 
one saddle point in the functional space. In this situation, the wave function which describes 
the phase slip process is determined so as to minimize the energy functional at each stage of 
the tunneling process. Of many tunneling paths, we adopt only paths which describe slips of 
2tt. Other paths of Air, 6ir, . . . can be omitted because of their much higher tunnel barriers. 

To determine the shape of the path, we temporarily fix the phase slip center between the 
ith and (i + l)th grain and restrict the average value of supercurrent, so add the following 
two terms in Hamiltonian Eq. @ 



where q is the "phase jump" at the phase slip center satisfying — tt < q < tt. The "winding 
number of the phase" n is a real number, but not necessarily an integer. The extremal 
solution of the energy functional with additional terms is 
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Hi = n, 9 i+1 -6i = k= _ 2 (ttJ), 

rij = n, 9j + i — 9j = q. (3) 

When we take the phase jump q as the collective coordinate describing each stage of the 
tunneling process, Eq. (Q) is the energy minimum configuration with given q and n. When 
q = k, Eq. (H) gives the local minimum solution {n m ^, m ,i}: 

^m,i = W>7 ^m,i+l 9 m i = "J^- (4) 

and when q = n — n, the saddle point solution {n S) i, 9 St i}: 

n 8>i = n, 9 Sii+1 - M = k= _^ - (i ^ j) , 

n St j = n, 9 Stj+1 - 9 sd = 7r - k. (5) 

Now the wave function is constructed which interpolate two local minimum points in the 
functional space. 

Hereafter we adopt k, rather than q, as the collective coordinate, because using k allows 
us to consider the whole of the tunneling process while q doesn't. The potential energy is 
now a function of the collective coordinate k, 

V (k) = -Ej [{N - 2) cos k + cos {2nn - k (N - 2)}] . (6) 

Next, we discuss the zero-point fluctuation of the collective coordinate k around the local 
minimum. As already mentioned, rij and 9i are canonically conjugate variables, therefore the 
field necessarily fluctuates around the local minimum configuration Eq. (|4]) and a certain 
mode of this fluctuations corresponds to the zero-point fluctuation of k which nucleates 
phase slips. In order to estimate frequencies of collective modes, we add the fluctuation 
terms {Srii, 89^} to the local minimum configuration Eq. (Q), substitute it into the Euclidean 
action, Eq. (|I|), and retain the terms up to the second order of Srii and 86 i, then we obtain 

(7) 




Here, G is the Green's function for 5ni and 59i, and Fourier component of G 1 is given by 
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where a is the distance between nearest neighbor grains. The poles of G give collective 
modes, whose dispersion relation is 



. i _ _ 9 ka 1 _ / 27m \ 9 fca 
hw k = 2 J {E G + E c sin 2 — \ E 3 cos — — sin 2 — . (9) 
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Fig. m suggests that a phase slip has a local character, thus as the effective frequency of the 
zero-point oscillation of k, we can extract from Eq. <M) the mode whose wave length is a: 



hu eff (K n ) = ftw k= i = 2y (E G + E c ) Ej cos (^^jj ■ (10) 

We have already determined the potential energy, Eq. (^), so the effective inertial mass 
of k can be estimated from the frequency of the collective mode and the curvature of the 
potential energy at the local minimum. In this procedure, we assume that the inertial mass 
is not greatly dependent on the collective coordinate, and approximate it to the value at 
local minimum. From Eq. (|6|) the potential curvature at the local minimum k ~ n n is 

V"(K n ) =E 3 (N- 2) (N - 1) cos (J^) , (11) 

where the prime denotes differentiation with respect to k. We can derive the effective inertial 
mass for k from Eq. ( |10D and Eq. (|TT|). 
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We now consider the tunneling problem from the local minimum K n to the next lower 
local minimum K n _i. According to WKB formula, the tunneling rate T(n n —>■ « n _i) is given 
by 

T(K n « n _i) = — — — exp <^ -~S(K n ) V , 



S{K n ) = J J J 2M(n n ) (v(k) - \hu k= ^j , (13) 



where k~ and k+ are the classical turning points, located between n n -i and K n . Here k+ 
doesn't coincide with n n because we adopt not the instanton method, but WKB method 
which explicitly includes the zero point energy in the action. This point will be discussed 
later. Besides, we implicitly assume the existence of the environments, so that the energy 
is conserved before and after the tunneling. 

In the limit of large N and small bias current, we can reduce the action in Eq. (|13D up 
to the first order of a small parameter to obtain 



S(K n ) ~ 2hJ x 



(14) 



where J = y E J^ Ec is the coupling parameter representing the strength of the phase coher- 
ence. F and E denotes the first and second kind of elliptic integrals respectively. 

From the tunneling rate, we can derive an expression for the resistivity. Combining 
Josephson's relation 

d6 tot 2e , 



dt h 

and the equation of motion of the phase 

dO tot 
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27iT(K n -> /« n _i), (16) 

we conclude that the voltage at metastable state k ~ K n is 

V{n n ) = — r(/t n -> K n -i), (17) 
e 

where 8 to t denotes the total difference of the phase throughout the array. The zero bias 
resistivity R is defined as follows: 

dV{K n 
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where I{n n ) is the applied current equal to E^sm.n n 



(18) 



In the limit N — > oo, we can actually calculate the zero bias resistivity in the analytical 
form. From Eq. ([13]), (^4|), ([XT]) and flT8|), i? is obtained in the following form: 

exp 

Note that J contains both Eq and Eq. Many theoretical works have treated only the limiting 
case in which either Eq or Eq is absent, while our derivation includes the case in which the 
two charging energies are comparable. Incidentally, it is reported that the resistivity is 
nonlinearly dependent on the applied current in the case of superconducting wires |§, but 
that might be caused by applying, instead of WKB method, the instanton method || which 
doesn't explicitly take into account the zero point energy in the action. This is the reason 
why we adopt WKB method, instead of the instanton method. 

The resistivity Eq. ([19]) is shown in Fig. |2] compared with the experimental result. 
Because the experiment is in the limiting case where Eq 3> Eq, the data are plotted against 
the reduced parameter Jo = J JK Without any fitting parameter, our theory and the 
experiment show good agreement asymptotically in the superconducting side (Jo > 0.5). 
However WKB approximation, on which our theory is based, is expected to be valid in the 
strong coupling limit (Jo ^> 1), so more experimental data in the stronger coupling regime 
(Jo > 0.8) are necessary to confirm our theory. Further, the data of arrays in which Eq and 
Eq are comparable are also hoped for. 

We can also estimate the S-I transition points for various number of grains N. In a strict 
sense, there is no critical point for a finite N because there is no sharp phase transition 
in a finite size system. But instead we can consider some critical region within which the 
crossover from the superconducting state to the insulating state takes place. In fact, in the 
experiment |D, the behaviors of resistance for finite iV's (63,127,255) change around Jo ~ 0.5, 
so we can state that the transitions occur around there. The positions of critical regions are 
slightly dependent on N. To be more specific, with increasing N, the crossover is seen to 
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occur at larger J. This seems to correspond to the ascent of the transition temperature due 
to the size effect in the case of a 2D XY model, which can be detected experimentally in 



superfluid He in mesoporous media [10, f f 



In our theoretical derivation, when the zero point fluctuation of k and it's potential 
barrier are comparable, the phase cannot be localized and the phase coherence is destroyed. 
So we can define this configuration as the boundary between superconducting and insulating 
states. The critical value Jq is estimated from Eq. @ and Eq. (|l0l), 

Jc ~ {(N - 1)-(JV -3) cos . (20) 



For large N, Eq. (pOj) reduces to 



2 \ 4 \nJ 

Eq. fl21~|) qualitatively explains the size dependence of the critical region in the experiment, 
although the quantitative agreement is rather poor. This quantitative discrepancy may result 
from the harmonic approximation we have employed for estimating zero-point fluctuation 
of K. 

Now, we give some comments to our argument up to this point. There might be an 
assertion that we should take into account the translational symmetry about the position of 
a phase slip center and therefore multiply the tunneling rate by the number of virtual centers 
of a phase slip formation, say, the number of junctions. But the mode we have adopted as the 
zero point fluctuation of k is created by the collective motion of the field and not localized at 
some junction, so the degree of freedom about the position of a phase slip center is already 
summed up. Besides, it is worth considering additional tunneling paths slightly shifted from 
the valley paths which we have adopted. This supplement will be crucial when the effective 
inertial mass of a phase slip is greatly dependent on the tunneling path and the correction 
by these paths probably increase the tunneling rate. On the contrary, it is widely known 
that the coupling of superconducting component with the environments generally suppresses 
the tunneling rate [|12|| , which we have ignored in this Letter. To evaluate the contribution 



from this dissipation effect is also left for the future work. 
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In conclusion, we have studied the quantum nucleation process of a phase slip in a ID JJ 
array using the collective coordinate method. The theoretical value of the zero bias resistiv- 
ity calculated by means of WKB method in the case of large N limit is in good agreement 
with the experimental data asymptotically in the superconducting side. We also have es- 
timated the S-I transition point for an array with a general number of junctions, which 
roughly agrees with the experiment. All of our results suggest that the intrinsic resistivity 
in the superconducting side in a ID JJ array at zero temperature can be understood with 
the concept of quantum nucleation of phase slips and it's value can be actually estimated. 

This work is supported by Grant-in-Aid for Scientific Research from the Ministry of 
Education, Science, Sports, and Culture of Japan. 
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FIGURES 




FIG. 1. The phase of the interpolation function Eq. (^) for the case ./V = 63, n = 2 and q = 3 
is plotted (solid line). The dotted line depicts the local minimum solution Eq. (||). 
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Jo = (£j/£ G ) 1/2 

FIG. 2. The zero bias resistivity plotted against the coupling parameter Jo. In the supercon- 
ducting side (Jo > 0.5), our theoretical value is in good agreement with the experimental data, 
especially with the ones for N = 127 and N = 255. 
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